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Abstract 

We present and study a family of metrics on the space of compact subsets 
of R N (that we call "shapes"). These metrics are "geometric", that is, they are 
independent of rotation and translation; and these metrics enjoy many interesting 
properties, as, for example, the existence of minimal geodesies. We view our space 
of shapes as a subset of Banach (or Hilbert) manifolds: so we can define a "tangent 
manifold" to shapes, and (in a very weak form) talk of a "Riemannian Geometry" 
of shapes. Some of the metrics that we propose are topologically equivalent to the 
Hausdorff metric; but at the same time, they are more "regular", since we can hope 
for a local uniqueness of minimal geodesies. 

We also study properties of the metrics obtained by isometrically identifying a 
generic metric space with a subset of a Banach space to obtain a rigidity result. 
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1 Introduction 

A wide interest for the study of shape spaces arose in recent years, in particular inside 
the Computer Vision community. 

There are two different (but interconnected) fields of applications for a good Shape 
Space in Computer Vision: 

Shape Optimization where we want to find the shape that best satisfies a design goal; 
a topic interest in Engineering at large; 

Shape Analysis where we study a family of Shapes for purposes of statistics, (auto- 
matic) cataloging, probabilistic modeling, among others, and possibly to create 
an a-priori model for a better Shape Optimization. 

To achieve the above, some structure is clearly needed on the Shape Space, so that 
our goals can be studied and the problem can be solved. 

Remark 1.1. Note that, for the purpose of Shape Optimization, shapes are usually 
intended "up to rotation, translation and scaling "; for this reason, when we wish to 
distinguish between the two, we will call a space for Shape Optimization a "preshape 
space ". 
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1.1 Shape spaces 

In general the "Shape Space" 2" will be a suitable choice of subsets of R N . 
A common way to model shapes is by representation/embedding: 

• we represent the shape A by a function ua 

• and then we embed this representation in a space E, so that we can operate on 
the shapes A by operating on the representations ua', 

for example, if E is a Banach space with norm || ■ ||, we can define a distance of shapes 
by d(A, B)=\\u A - u B \\. 

Most often, this representation/embedding scheme does not directly provide a Shape 
Space satisfying all desired properties. In particular, in many cases it happens that the 
representation is "redundant", that is, the same shape has many different possible rep- 
resentations. An appropriate quotient is then introduced. 

There are many examples of the representation/embedding/quotient scheme in the 
literature; for the case of generic subsets of R N , 

• a standard representation is obtained by associating a closed subset A to the 
distance function 

ua{x)= inf \x - y\ (1) 

y€A 

or the signed distance function 

b A (x)=u A (x) - u r n\ A (x) (2) 

We can then define a topology of shapes by deciding that A n — > A when ua„ — > 
ua uniformly on compact sets. This convergence coincides with the Kuratowski 
topology of closed sets. 

We can also operate "linearly" on shapes by operating on ua or 6^: so we can 
define shape averages and shape principal component analysis. Note that in 
general a linear combination of (signed) distance functions will not be a (signed) 
distance function: so any linear operation must be followed by an ad hoc correc- 
tion. For example, given two shapes Aq, A%, we can define an interpolation A t 
for t e [0, 1] by setting A t = {x \ tb Al (x) + (1 - t)b Ao (x) < 0}. 

This Shape Space is not independent of the position: when it is used for shape 
analysis, a registration of the shapes to a common pose is often performed (but, 
see also sec. 12.1. 11 1. 

• A. Duci et al (see 17] 0) represent a closed planar contour as the zero level of a 
harmonic function. This novel representation for contours is explicitly designed 
to possess a linear structure, which greatly simplifies linear operations such as 
averaging, principal component analysis or differentiation in the space of shapes. 

• Trouve-Younes et al (see [9 1, [27 1 and references therein) modeled the motion 
of shapes by studying a left invariant Riemannian metric on the diffeomorphisms 
of the space R"; to recover a true metric of shapes, a quotient is then added. 

But the representation/embedding/quotient scheme is also found when dealing with 
spaces of curves: 
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• In the work of Mio, Srivastava et al. l20l |T9j HU, smooth planar closed curves 
c : S 1 — > R 2 of length 2ir are parametrized by arclength and represented by the 
angle function a[0, 2tt] — > IR such that 

c(s) = (cos(a(s)), sin(a(s))) 

then the angle function is embedded in a suitable subspace N of L 2 (0, 2tt) or 
W 1,2 (0, 2tt). Since the goal is to obtain a Shape Space representation for Shape 
Analysis purposes, a quotient is then introduced on N. 

• Another representation of planar curves for Shape Analysis is found in Younes 
||3TI . In this case the angle function is considered mod(7r). This representation 
is both simple and very powerful at the same time. Indeed, it is possible to prove 
that geodesies do exist and to explicitly show examples of geodesies. 

• Metrics of "geometric" curves (that is, curves up to the choice of parametriza- 
tion) have been studied by Michor-Mumford [ 1 8] [TT] [16] and Yezzi-Mennucci 
Ej] |28] [29); more recently, Yezzi-Mennucci-Sundaramoorfhi ll26l l25l l24l |231 
[T5l |2T1 l22l have studied Sobolev-like metrics of curves and shown many good 
properties for applications to Shape Optimization; similar results have also been 
shown independently by Charpiat et al [5 1. 

Remark 1.2. In this case, shapes are modeled as immersed parametric curves 
c : S 1 — > R N , for the sake of mathematical analysis; a quotient w.r.t the group 
of possible reparametrizations of the curve c (that coincides with the group of 
diffeomorphisms Dif^S* 1 )) is applied afterward to all the mathematical struc- 
tures that are defined (such as the manifold of curves, the Riemannian metric, 
the induced distance, etc.). 

1.2 Goals 

We remarked that, in Shape Analysis, shapes are usually considered "up to rotation, 
translation and scaling", but even in Shape Optimization, to a certain degree, our theory 
should be independent of rotation and translation: that is, whatever we do with shapes 
should not depend on "where in the plane" we do it. 

In the rest of the paper we will denote by X the family of the nonempty compact sets 
in R N , and we will build many examples of metrics d on I. We will always require 
these metrics to be euclidean invariant. If A is an euclidean transformation of the 
space (a rigid transformation), then 

d{Ml x ,AD. 2 ) = d(Oi,Q 2 ) • 
What other properties may be interesting for applications? 

As mentioned before, a goal of Shape Optimization is to define shape metrics, 
shape averages, shape principal component analysis, shape probabilities. . . 

For example, if we represent shapes Aj,j — 1 ... n by their signed distance func- 
tion bAj , then we may define Signed Distance Level Set Averaging 

1 N 

A={x\f(x) = 0j, where f(x) = ^J^-CO (3) 

n=l 

A benefit of this definition is that it is easily computable; a defect is that, if the shapes 
are far way, then A will be empty. Another defect is that this definition is quite ad hoc: 
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it is not coupled with any other structure that we may wish to add to the Shape Space, 
such as a metric d. We may then look at the problem in the other direction. 

Considering a generic metric space (M, d), define the Distance Based Averaging^ 
of any given collection <ii . . . a n € M, as a minimum point a of the sum of its squared 
distances: 

n 

a = arg min \_. d(a, aj) 2 (4) 

a . 

i=l 

Supposing now that the Shape Space X is given a metric d, we can use the abstract 
definition above to define shape averages; this definition has many advantages. Namely 

• it comes from a minimality criterion, so it is "optimal" in a certain sense (con- 
trary to the definition (f3]l). 

• If the distance is invariant w.r.t. a group action, then the shape average is as well 
(see sec. l2.1.lT ). For example, in the case of geometric curves, where the distance 
is independent of parametrization, then the shape average will be independent 
of the parametrization of a\ . . . a n . 

• It coincides with the arithmetic mean in Euclidean spaces; more in general, when 
X is a smooth submanifold of a Banach space and a\ . . . a n are near enough, then 
a is an approximation of the arithmetic mean. 

In particular, the average of two shapes A\, A2 is the midpoint, that is a shape A 
such that 

d(A 1 ,A)=d(A 2 ,A) = ^d(A u A 2 ) 

We are then, however, bound by this result (whose base idea goes back to Menger - 
see sec. §4.i. 1 in JT4] for more details) 

Theorem 1.3. Supposing that d is complete and intrinsic (see sec. \2.1h then the fol- 
lowing facts are equivalent: 

• for any two shapes A\ , A 2 there is a midpoint; 

• for any two shapes A\ , A 2 there is a minimal geodesic connecting them. 

For this reason, we end up studying whether the Shape Space admits minimal 
geodesies (in theorem r3.18b . 

1.2.1 Tradeoff 

Unfortunately a tradeoff (that is well known in the Calculus of Variations) arises; 

• on one hand, a Shape Space that is useful for Shape Optimization should possibly 
be equipped with a topology that makes functionals "regular", so that suitable 
minimization methods can be used; to this end, the topology should have many 
open sets. 

• On the other hand, to prove existence of average points and of geodesies (that 
are useful in Shape analysis), it is sufficient that certain bounded sets be compact 
(cf. 12.21 1231 12.4b : to this end, the topology should have few open sets. 

'also known as Karcher mean, but it is also sometimes attributed to Frechet, in 1948 
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We can exemplify this as follows. 

• As we mentioned before, it was shown in Il25ll24l that flows for Active Contour 
methods that use a Sobolev metric are more robust to noise and converge faster 
than standard flows. To explain the rationale, suppose H, L are two metrics, H 
being stronger than L. When evolving the shape A, the L-flow will move to- 




wards a shape B with small scale deformations (such as those induced by noise), 
since B is nearer to A in the the i-induced distance; whereas the ii-flow will 
move towards the shape C with large scale deformations, since C is nearer to A 
w.r.t its related distance. 

• Suppose now though that a dataset contains a template shape A; an algorithm is 
given a version B of A that was corrupted by noise, and different shape C, and 
it has to decide what is the best match to A. In this case, the weaker metric L 
would associate A to the correct shape B, whereas an algorithm using the metric 
H would fail to associate A to B. 

For all those reasons, it is quite difficult to find a Shape Space that is suited both 
for Shape Analysis and for Optimization. 

1.3 Plan of the paper 

The plan of the paper is as follows: we foremost provide base definitions, and we 
propose some results in the theory of metric spaces, in particular when they are iso- 
metrically embedded in Banach spaces. Considering the space I of compact sets we 
review the well-known Hausdorff distance, and its properties; we successively propose 
a class of metric spaces that are similar to the Hausdorff distance, while at the same 
time enjoying some extra properties that may be useful in applications. 
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2 Metric spaces and embeddings in Banach spaces 



2.1 Metric spaces 

We recall some basilar definitions and results in the abstract theory of metric spaces. 
Suppose that (M, d) is a metric space. We induce from d the length len d 7 of a contin- 
uous curve 7 : [a, 0] — > M, by using the total variation 

n 

\en d 7^ sup V dfrfa-i), l(U)) , (5) 

where the sup is carried out over all finite subsets T = {to, • • • , i„} of [a, f3] and 
to < • • • < t n . 

We define the induced geodesic distance d 9 by 

d 9 {x,y)=mi\en dl j, (6) 

7 

where the inf is taken in the class of all continuous curves 7 connecting x to y. If the 
inf is a minimum, the curve providing the minimum is called a geodesic. Note that it 
may be the case that d 9 (x, y) = 00 for some choices of x, y. Note also that d 9 > d. 

We may consider the metric space (M, d 9 ); but the topology of (M, d) and (M, d 9 ) 
may be quite different, as we see in this example: 

Example 2.1. Let us consider the subset M=tJj(E) ofR 2 , where 
£= f [0,l] x ({0}U{l/n|n>0}) 

and 

( See fig. [7]). We associate to M the distance d induced by the usual distance of R 2 
onto M and d 9 for the geodesic distance. Then we have that (M, d) is obviously 
compact whereas [M, d 9 ) is not: indeed X n =ip(\, 1/n) does not admit a converging 
subsequence, since d 9 (x n , y m ) = 2 for all n 7^ m. 



Figure 1: example [2TTI 

When d = d 9 , we will say that the metric space is path-metric, or that d is intrin- 
sic. 

Note that the length len d defined by d 9 coincides with lcn d , and then d 9 = (d 9 ) 9 : 
d 9 is always intrinsic. 

We will use the following proposition: 

Proposition 2.2. if for a choice of p > 

W{x,p)Z{x\eP(x,y)<p} (7) 

is compact in the [M, d) topology, then x and any y € D 9 (x, p) may be connected by 
a geodesic. 
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The proof is simply obtained by the direct method in the Calculus of Variations (see 
Thm. 4.24 in |fl4l ). 

We also state these simple propositions. 

Proposition 2.3. If for ax G M for all choices of p > 

®(x,p)={x | d(x,y) < p} 
is compact then (M, d) is complete. 

Proposition 2.4. Suppose that a\ . . . a n € M are given; a sufficient condition for the 
existence of the Geodesic Distance Based Averaging aof a\ . . . a n 

n 

a = argmin Q r(o) , where r(a)= / ^ d 9 {a, a 3 ) 2 (8) 

.7=1 



is that, defining 



p* = min r{ai) 

i=l,...n 



we have that p* < oo and that D 9 (ai, 2y/p* + s) is compact in the (M, d) topology, 
for s > small. 

Proof. Note first that the infimum of r(a) is finite, since it does not exceed p* . Recall 
that 

d 9 (a, aj) = inf L 

where lj is the length of a Lipschitz curve 7$ connecting a, ai. So we can rewrite the 
problem © as 

n 

inf 0(71... 7„), where 0(71 .. -7n)=y^Gi) 2 

71 ■■■In * ' 

3=1 

where the infimum is computed on all choices of Lipschitz curves 71 ... j n of length 
1 1 . . .l n connecting Oj to a common point x E M; for simplicity we represent them as 
7i : [0, Z 2 ] — > M parametrized by arc parameter. By triangular inequality 

d 9 (a ul3 {t)) < d 9 (a l ,x)+d 9 {x ll] (t)) < k + lj 

Let then 7^ be a sequence of choices that converges to the infimum: 

0(7i,fe ■ ■ ■ 7n,fc) -^k inf 0(71 . . . 7„) 

71 —7n 

so for large fc, 

0(7i,fc ■ ■ -7n,fc) < P* + e 



but then in particular ^ < ^/p* + e hence 



d s (ai, 7j (t)) < 2^/^+1 

so all the curves are contained in a compact set. By Ascoli-Arzela theorem, we can 
then extract a uniformly convergent subsequence, and use the fact that the length is 
lower semi continuous. □ 
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A similarly proposition can be stated for d: 
Proposition 2.5. Suppose that a%. . .a n S M are given; let 

n 

p* = min y d(a,i, a,-) 2 (9) 
i * — ' 
3=1 

and i* the index that achieves the above minimum: suppose that D(<2j*, ^Jp* + e) 
is compact for e > small: then there exists a point a that is the Distance Based 
Averaging of at... a n , as defined in @. 

2.1.1 Distances, quotients and groups 

Let dAi(x, y) be a distance on a space M, and G a group acting on M; a distance ds 
may be defined on B — M /G by 

d B ( [x] ,[y])= inf d M (x,y)= inf d M (gx, hy) 

xe[x],ye[y] g,h<£G 

that is the lowest distance between two orbits; we write ds (x, y) for simplicity. 
If d^j is invariant w.r.t. G, i.e. 

dM{gx,gy) = d M (x,y) Vg E G 

then 

d B {x 7 y) = inf d M {gx 7 y) (10) 

It is easy to see that ds satisfies the triangular inequality; but it may be the case that 
ds(x, y) = even when x ^ y. We state a simple sufficient condition 

Lemma 2.6. If the orbits are compact, then ds is a distance. 

When studying metrics dona Shape Space X, the quotient is particularly useful in 
at least two cases: 

• when we want to pass from a preshape spaced to a shape space: in this case, G 
is the Euclidean group of rotations and translation (and sometimes of scaling); 

• when the representation is redundant: for example, in remark 1X721 we would set 
G = Dif^S" 1 ). 

2.2 Embeddings in Banach spaces 

In most of what follows, we will be able to identify M (using an isometry i) with 
a subset of a Banach space E. We remark that an isometry is a map i such that 
d(x,y) = \\i(x)—i(y)\\ (and this should not be confused with the concept of isometrical 
embedding of Riemannian manifolds). 

2 cf. remark fOI 
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2.2.1 Radon-Nikodym property 



The following result from [ 1 ] will come handy: 

Theorem 2.7. Suppose that E is the dual of a separable Banach space. Letj : [a, b] — > 
E be a Lipschitz curve; then, by thm. 8.1 in SJ]/, for almost all t there exists the deriva- 
tive j(t) € B that is defined as 

7(f)^-lim 7(t + T) ~ 7(t) (11) 
where the limit is done according to the weak-* topology; and moreover, 

7(* + t)-7(*) 



117(011 = 



(12) 



so ||7(i)|| coincides with the metric derivative, that is studied in ft2§. 
There follows easily (by applying scalar products to ( 11 11 1) that 

7(6)- 7 (a) = f j(t)dt (13) 

J a 

and 

,6 

len d 7 = / ||7(<)ll* (14) 

J a 

( this last by thm. 4.1.1 in /E]?). 

It is common to say that E enjoys the Radon-Nikodym Property, when the limit 
in (fTTT i exists in the strong sense, and for almost all t. Note that the Radon-Nikodym 
Property does not hold in general: consider the map t i— ► l^t+i] in L 1 (R), whose 
derivative should be t i— ► <5t+i — <5f. 

We now recall this basilar definition: 
Definition 2.8. a Banach space E is uniformly convex if Ve > 035 > 0, 

Vx,yeE,\\x\\ < l,\\y\\ < ljz-vll > e => || + l/)/2|| < (1 - 5) . 

Examples of uniformly convex Banach spaces include L p (fl, A, fi) forp 6 (1, oo). 
Uniformly convex Banach spaces have many interesting properties: for example, they 
are reflexive (Milman Theorem, III. 29 in [4|); moreover, if x n — > x in weak sense and 
limsup ||x„|| < |jx|| then x n — > x in the strong sense (prp III. 30 in 101). 

So we obtain a sufficient condition: 

Corollary 2.9. ;/ E is uniformly convex and separable, then it enjoys the Radon- 
Nikodym Property f indeed eqn. (11 It and eqn. (1121 l imply that the limit in (11 11 1 is va/iof 
also in the strong sense). 



2.2.2 Embeddings in uniformly convex Banach spaces 

If E is uniformly convex then in particular the closed ball {x | ||x|| < 1} is strictly 
convex; this has a curious implication. 

Lemma 2.10. Suppose the closed balls in E are strictly convex. Consider E as a 
metric space, with distance ds(x, y) = \\x — y\\. The segment connecting x,y G E is 
the unique minimal geodesic (up to reparametrization). 
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Proof. We will prove that, for x, y, for any minimal geodesic 7 : [0, 1] — > M connect- 
ing x to y, if 7 is reparametrized to arc parameter then 7(1/2) = (x + y)/2; iterating 
this reasoning with finer subdivision we obtain that j(t) = (tx + (1 — t)y). 

With no loss of generality, up to translation and scaling, suppose y = —x and 
|| a; || = 1. The segment t 1— > tx is a geodesic for t 6 [—1, 1], by the theorem [2.71 
and its length is 2. Suppose now that 7 : [—1,1] — > M is another geodesic: then 
len7 = 2, and, up to reparametrization, ||7| = 1 at almost all points; in particular, 
setting z = 7(0), \\z — y\\ < 1 and ||x — z\\ < 1; but then, by triangular inequality, 
\\z + x\\ = \\x — z\\ = 1. Suppose that z ^ 0; then || (z + x) — (x — z)|| > 0; by strict 
convexity, though, this implies that ||((z + x) + (x — z))/2\\ = ||x|| < 1, and this is a 
contradiction. □ 

Theorem 2.11. Suppose that (M, d) is a complete space, and that i : M — > E is an 
isometrical immersion in a uniformly convex Banach space E. If given x, y G M, 
d(x, y) — d 9 (x, y), then the segment connecting i(x), i{y) is all contained in i(M). 

In particular, if (M, d) is path-metric then i(M) is convex, and then any two points 
in M can be joined by a unique minimal geodesic ( unique up to reparametrization). 

Proof. Note that i(M) is complete, and then it is closed in E. We will prove that, for 
any x, y € i(M), [x + y)/2 6 i(M); we can then iterate this idea to further subdivide, 
and since i(M) is closed then this proves the whole segment connecting x, y is in i(M); 
for the above lemma, the segment is the unique minimal geodesic. 

We now fix x, y £ i(M): there must be paths j n : [—1, 1] — > i{M) connecting x 
to y with length len(7) < L„=||x — y\\ + 2/n. 

As in the lemma before, we suppose for simplicity that y = — x and ||x|| = 1 
(so L n — 2 + 2/n); and we reparametrize so that ||7 n || = 1 + 1/n: hence setting 
z n = 7„(0) 

\\z„ +x\\ < 1 + 1/n , \\x - z n \\ < 1 + 1/n . 
and then by triangle inequality \\z n + x\\ —> l,\\z n — x\\ —* 1. Setting 

W n = (Zn +x)/\\Zn + x\\ , V n = (x - Z n )/\\z n - x\\ 

we can prove that || (w n + v n )/2\\ — > 1 hence by the uniform convexity of E we obtain 
that w n — v n — > and then z n —* 0. Since z n 6 i(M) and i{M) is closed then 
e i(M). □ 

The above is a "rigidity theorem", in that it restricts the class of metric spaces that 
can be isometrically embedded in a uniformly convex Banach space E. 

Corollary 2.12. a complete compact finite dimensional Riemannian manifold M can- 
not be isometrically embedded in a uniformly convex Banach space E: indeed in this 
space M there are two points that can be joined by more than one minimal geodesic. 

When E is not uniformly convex, on the other hand, strange behaviours arise. 

Remark 2.13. Let L°° = L°° (p,, A, fi) and suppose 57 is not an atom of /i, that is, 
suppose the dimension of L°° is greater than 1. Given generically /, g E L°°, there is 
an uncountable number of minimal geodesies connecting them. 

Proof. We can assume without loss of generality that g = and that ||/|| = 1, Let 
A = {|/| = 1}. We will prove that if there is only one geodesic then |/| = 1^. 
Indeed if |/| ^ t A then < 1} > 0. Let < t < 1 be such that fi{\f \ < t} > 0; 
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obviously fi{\f \ > t} > since ||/|| = 1; let A' = {\f \ > t} and A" = {\f \ < t}. 
Given any diffeomorphism b : [0, 1] — * [0, 1] with b'(s) < l/t, 

j(tptfl A , + b(t)fl A „ 
is a geodesic. Indeed its derivative is 

7 '(*)=/lA' + b\t)ft A „ 

and || 7' (t) || = 1 by construction. 

The family of / s.t. |/| = A1a is closed and has empty interior. □ 

The idea of isometrical embedding is quite powerful: indeed any separable metric 
space may be isometrically embedded in l°° (that is the dual of the separable space 
i 1 ): so the breadth of application of the theorem [277] is general, and is at the basis of 
many results in [ 1 1. But the embedding in £°° that is studied in [ 1 1 is not suited for our 
practical applications: 

• it would not respect the geometric properties of the space (as we discussed in 
sec. II. 2\ 

• it would be too difficult to find a satisfactory notion of "shooting of minimal 
geodesies" using this embedding. 

For all above reasons, we will consider isometrical embeddings in this paper as well 
but we will (for the most interesting applications) use an explicitly chosen embedding 
in uniformly convex Banach spaces. 

2.3 Definitions 

We introduce some definitions that will be used in the rest of the paper 
We will write sVt = max{s, t} and s + — max{s, 0}, when s,t€ R. 
We will write B(x, r) or B r (x) for the open ball of center x and radius r > in 

R N ; we will shortly write B r for B r (0). Similarly D r (x) will be the closed ball of 

center x and radius r > in R N , and D r — D r (0). 
We define the fattened set to be 

A + D r = {x + y I x G A, \y\ < r} = (J D r (x) = {y \ u A (y) < r}. 

xeA 

This fattened set is always closed, (since the distance function u A {x), that was defined 
in (Q~J|, is continuous). 

We will say that a family A.- Le i of sets in R N is equibounded if there is a R > 
such that A.- L C Dr for all i. 

We denote by C N the N dimensional Lebesgue measure, and ^^=£^(£1); we 
write shortly J A f(x)dx for the Lebesgue integral. 

2.4 Hausdorff distance 

A fundamental example of metric on X is the Hausdorff distance 

d H (n, il')= inf{<5 >0|ft'c (n + D s ),nc (0' + D 5 )} . 
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It is not difficult to verify that 

d H (Q,tt')= sup \uq(x) - uq'(x)\ , (15) 

x£R N 

see for example Thm. 2.2 in ch. 4 in Delfour-Zolesio |6|. 
This metric enjoys many important properties. 

Theorem 2.14. The metric space (T 1 du) satisfies: 

• given r > 0, the family of r-bounded compact sets 

{fl e 2 | fl C A} 
is compact; in particular, the set 

D={n | d H (n,n r ) < P } 

is compact; 

• (I, du) is path-metric (that is, dn — {dn) 9 ) 

• consequently, by Prp. \2.2\ any two £l,Q' 6 X may be joined by a minimal 
geodesies; 

• and moreover, by Prp. \2.3\ (X, dn) is complete. 

The first statement is a well known property of the Hausdorff distance, see e.g. 
|6| pag. 194. By exploiting the characterization ( fT3T l. it also follows from a diago- 
nal/compactness argument and the following rigidity property: 

Lemma 2.15. Let Q n be closed sets, and suppose that lim n uq u (x) = f(x) for all x 
in a dense subset D ofR N . Then there is a closed set fl such that Ufi(x) = f(x) for 
all x € D, and u$i„ — > u$i uniformly on compact sets; moreover if (and only if) f2„ is 
equibounded then UQ n — > un uniformly. 

Proof. The proof may follow from the theory of Viscosity Solutions: it is well known, 
indeed, that uq is the unique solution to a properly defined Eikonal equation; and that 
viscosity solutions do enjoy the required rigidity property. 

We propose here instead a direct proof. We set u n =Un n ; it is easily proved that u n 
is 1-Lipschitz, that is 

\u n (x) - u n (y)\ < \x - y\ Vx,y (16) 
so passing to the limit in the above ( TTST ). we obtain 

\f(x)-f(y)\<\x-y\ Vx,yeD (17) 

and then there is an unique extension of / to a positive function g : R N — > IR that is 
again 1-Lipschitz, that is, 

\g(x)-g(y)\ <\x-y\ Vx,y . (18) 

It is easy to prove that u n (x) — » g(x) for all x, and actually (by imitating the proof 
of Ascoli-Arzela theorem) that u n — > g uniformly on compact sets. 
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Let fi = {g = 0}; to conclude the proof, we need to prove that g = uq. To this 
end, we first prove that g > uq: indeed, fixing x, u n (x) = \x — y n \ for at least one 
point y n G fi n ; since u n (x) — ► g(x), then the sequence {y n } is bounded, so (up to 
a subsequence rife) it converges to a point y; since the family u n is 1-Lipschitz and 
u n (y n ) = then g(y) = 0, that is y G ft: hence 

g(x) = limu„ fe (x) = lim \y nk - x\ = \y - x\ > u n (x) . 

k k 

Conversely, let y G ft be such that uji(-t) = \x — y\; then by (TT~8T > g(x) < g(y) + \x — 

y\ = \ x - v\ = un(x). 

To conclude, supposing that fl n is equibounded, then choosing R > such that 
fl n C Da, we know that un n — ► Ufj uniformly on D^, so given e > for n large 
|u„ — u| < e in Z3/j and then 

u n (x) = inf (|x - y\ + u n (y)) < inf (\x - y| + + e) = u(x) + e 
y£D„ yeD R 

where the first and last equalities are due to the Dynamical Programming principle; 
and similarly we obtain that u(x) < u n (x) + e. The "only if" part follows from 
O- □ 

To prove the above second property in 12. 141 we may use the first property and the 
following Menger convexity result 

Proposition 2.16. Let A, B € T be two compact sets, then for all A € [0, l]there exists 
a compact set C such that 

• d H {A,C) = \d H {A,B), 

• d H (B,C) = (l-X)d H (A,B). 

Proof. We write [A] r = A + D r for the fattened set. Let /i = da (A, B). We consider 
the set 

C= {z\3x eA iV eB, \x -y\<fi, \x -z\< X/j,, \y - z\ < (1 - A)^} . 

We prove that the set C has the properties we need. In particular it is enough to prove 
only the first one because of the symmetry in the two conditions. If x G A then there 
exists y G B such that \x — y\ < p., then the a point z = (1 — \)x + Xy satisfies 

\x - z\ < Xp, \y-z\<(l- X)/d. 

Such a z must be an element of C and so we found an element of C with distance less 
or equal to A/i. This means that x G [C]a m and it is true for all x G A so A C [CJa^. 

Let's take now z G C. From the definition of the elements of C we have that there 
must exists x G A and y G B such that \x — z\ < Xfi, \y — z\ < (1 — X)p. This means 
that z G [A] Xf j,. This is true for all z G C so C G [A] Xll . 

To finish the proof we have to show that the set C is compact. It is clearly bounded 
because it is contained by We have to show that it is closed. Suppose we have 

a sequence {zfcjfe C C such that Zk — > z. Then for each Zk we can find two elements 
Xk G A, yk G B with the properties: 

\xk ~ Zfc| < Xp, \y k - z k \ < (1 - X)[i. 
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The sets A and B are compacts so we can chose a subsequence (for simplicity we use 
the same index k) such that Xk —> x E A and y^ — > y E B. It is obvious to see that the 
points x, y, z satisfy the following inequalities: 



\x - y\ < j«, \x- z\< A/U, \y - z\ < (1 - A)/x. 
This means that z is an element of C and this concludes the proof. 



□ 



Unfortunately (X, djj) is quite "unsmooth", as shown by this example (that is sim- 
ilar to !2. 131 — and for a reason!). 

Example 2.17. There are choices of el that may be joined by an uncountable 
number of geodesies. 

In fact we can consider this simple example: 



A = {x = 0, < y < 2} 

B = {x = 2,0 < y < 1} 

C t = {x = l,0<y< §}U 
{tf = 0,l<x<t} 

with 1< t < \/5/2: 



A 



Ct 

i 

I 



B 



and in the picture we represent (dashed) the fattened sets A + B^, 2 and B + B^, 2 . 

Note that d H (A,B) = V5 while d ff (A,C t ) = d H {B,C t ) = y/S/2: so C t are all 
midpoints that are on different geodesies between A and B. 

We conclude with a family of nice properties. 

Proposition 2.18. 1. The fattening map A i— > A + D\ is Lipschitz (of constant 
one). 

2. Given A > 0, the "fattened area map" L\(A)=C N (A + D\) is continuous. 

3. Consequently, the area map L(A)=C N (A) is upper semi continuous. 

4. Let # : X — > NUoobe the number #17 of connected components of a closed 
set f2. Then # is lower semi continuous in the metric space (X, dn)- 

As a corollary, the family of connected compact sets is a closed family in (X, dn)- 
Proof. 1. Obvious. 

2. if A„ — > A then for fixed e > and definitively in n, 

A n cA + D E , AcA n + D e 

and then 

An + D x C A + D £+x , A + D x -e C A„ + D\ 
passing to Lebesgue measures, 

C N (A+D X - e ) < \immiC N {A n +Dx) < \imsup C N (A n +D x ) < C N {A+D £+X ) 

" n 

and we let e — > 0. 
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3. Since it is the pointwise limit L\{A) [ L(A) for A — * 0. 

4. See Thm. 2.3 in ch. 4 in (6). 



□ 



3 L p -\ike metrics of shapes 

The definition of the Hausdorff distance by eqn. (15[ leads us back to the paradigm 
of representation/embedding; but in this case it is unfortunately not precise, since the 
Banach metric that we use, namely 

11/11 = 11/1100=^1/(1)1 

X 

is usually associated to the spaces Cb(R N ) of bounded functions — whereas the dis- 
tance function ua is not bounded! What follows is a simple yet effective workaround. 

Definition 3.1. We fix p G [1, oo]; we fix a function tp : [0, oo) — » (0, oo) monotonically 
decreasing and of class C , such that 

<p(\x\) € W(R N ). (19) 
Note that, for p < oo, the above is equivalent to asking that 

/■OO 

/ t N - l <p{t) p dt < oo (20) 
Jo 

and it implies that lirut_>oc tp(t) — 0;forp — oo we instead ask that lim^oo tp(t) = 
as an extra hypothesis. 

An example of such a function is <p(t) = exp(— t), or tp = (1 + t)~( N+1 " p . 
We will often write va — tp o ua for simplicity. 

Lemma 3.2. Let O C R N be closed and non empty; suppose p < oo; then 

(a) v n G U>(R N ) if and only if 

(b) f2 is bounded (and then f2 is compact). 

Proof. We first prove that (a) =£■ (6) by contradiction. Let us assume that 57 
is unbounded. Then there exists a sequence {x^} C £1 such that |at^| — > oo and 
d(xk, x q ) > 2 for all k, q G N, fc 7^ q. The sequence of sets Bi(xk) is disjoint. It is 
easy to see that Vn(x) > ^(1) for a; G (J fe Bi(xk), and then u ^ L p . 

Then we prove that (6) => (a). If is bounded we can find a ball Br such that 
fi C Br. Then easily we have uq > ur => Vq < vr, but »g £ L p (as is easily 
proved by Vr(x) = <p((\x\ — R) + ) and by (1201 )) and then also vq 6 LP . □ 

Definition 3.3. Given A, B G I, we define 

dp,tp(A, B)=\\tp(uA) - <p(ub)\\lp{rX) 
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By the above lemma, this distance is finite. We will often write d for d PtV in the 
following, for simplicity. 

The above distance is obtained by the representation of a shape A as Va, combined 
with the embedding of va in L P (R N ). For this reason, we may identify our shape space 
with 

n c ={v q \nei} (2i) 

that is a subspace of LP . 

Remark 3.4. By the definition of d, the map £1 1— > vq is an isometrical embedding of X 
inside LP, and the image is N c ; N c is a closed subset of LP, by the completeness result 
13.1 ll that we will prove in the following. 

We will exploit this embedding in the following, as in 33.61 
It is immediate to verify that d Pjtp satisfies these properties. 

• The embedding A i— > va is injective: if va=vb then ua=ub (since ip is mono- 
tonically decreasing, and so it is injective); but, by lemma l2~T5l this implies that 
u a = ub and then A = B; consequently, for all A, B £ 1, d PtV (A, B) = iff 
A = B. 



• d PiV is euclidean invariant, as we requested in sec. 11.21 

d P)¥! (fix, fi 2 ) < llufijz* + ||«n a \\lp- (22) 

forp < co, and 

(n 1 ,n 2 )<¥»(Q). 

Proof. When p < oo, by the Minkowski inequality we have that d(Qi, f2 2 ) < 
II vq 1 || lp + || vq 2 || lp ; moreover equality would hold only if vq 1 = —vq 2 and this 
is impossible; when p — oo we use the fact that ip > 0. □ 

• (Separation at infinity) given two bounded sets f2 2 we have 

lim dp !V> (n 1 ,rt 2 +T) ^\\vnA\ LP + \\vq 2 \\ L p ; (23) 

|t| — *-oo 



forp < oo, and 



lim d^ni, n 2 + t) = tp(0) ; (24) 

|r | — >oo 



Proof. For the case p < oo this comes from a general result for LP functions; 
for p = oo it derives from the hypothesis limt_ s . 00 ip(t) = 0. □ 

• (Scaling) If p < oo and A > is a rescaling of the space, then the rescaled 
distance may be expressed as 

^(Afii,An a ) = x N/p d p ^(n 1 ,n 2 ) (25) 
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where (p(r) = ip(Xr); indeed 



dp,ip(Xflx, AS1 2 ) P = / \v\nA x ) - v X n 2 (x)\ p dx (26) 



= X N J \vxn t (Xz) - v Mh (Xz) \ p d z (27) 

= X N J \<p(Xu ni (z))-<p(Xun 2 (z))\ P dz (28) 
= X N d p ^(n 1 ,Q 2 ) p (29) 




where to go from (|26"T i to (l27l i we used the change of variable x = Xz and the 
property of the distance function 



to change (|27j» to d28i 

Remark 3.5. The inequality (l22l i easily implies that the balls of the distance d in general 
are not compact sets. Indeed it is enough to consider a set 51 and the following ball: 
D = {A | d(A, 51) < 2r} with r = ||«f2||z,p. Then the sequence: {51 + nr} ngN with 
t e R N \ {0} is contained in D and it does not have any convergent subsequence. 

To continue with our study of d, we prove this fundamental inequality. 

Lemma 3.6 (Local equiboundedness). There is a continuous and increasing function 
b : R + —> R + with 6(0) = such that, for any Cl, Cl' £ X satisfying 



To prove the proposition for p € [l,oo), suppose that xo G 51', but xq K; for 
y G B(xq, r/2) recall the simple triangular inequality 



u\n(Xz) = Xun(z) 



(30) 



\\vq - vw\\lp < b(r), 



then 51' C 51 + D r . 



Proof. Set K=Q + D r . It is easy to check that 



v K (x) = tp((un(x) - r) + ). 



ua(v) > r - \x - y\ > \x - y\> u a >{y) 



hence 



vn(y) < f{r - \x - y\) < p(\x - y\) < vn>(y) 





\<P(\%0 - V\) ~ ¥( r - \ x o - V\)\ p dx = b(r) 



where 




r/2 
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where lun is the volume of the ball B±. It is easy to prove that b is continuous and 
increasing (by direct derivation); that 6(0) = and that lim r _ >00 b(r) = \\lp{\x\)\\lp- 
The case p = oo is simpler: in this case we can note that 

\\vn - ffi'lU > Vfi'{xo) - vn(x a ) > ip(0) - ip(r) 

and set b(r) — <p(0) — <p(r). □ 

Corollary 3.7. As a corollary we obtain that for d(il, CI') small enough, 



Remark 3.8. The above does not hold for arbitrarily large distance d(Cl, fl'): indeed, 
let O = {0} and 0„ = {nei}: then e?(Sl,i7„) — > 2||<^(|x|)||lp (as we mentioned in 
423}). 

We can also obtain a converse inequality, as follows 

Lemma 3.9. There is a family of continuous functions f R : [0, 1] — > R + with f R (0) = 
0, such that, for any fi, f2' €E X, //T2 diameter R and dfr(f2, < 1, then 



Proof. We provide the proof forp < oo. Note that if has diameter i? and dn < 
1, then il' has diameter at most R + 2. Up to translation, suppose that £?2i?+4 contains 
both Q and f2': then 

wn(a;),wn'(x) < <p(Qx\ - 2R- 4)+) 

so 

/ Mx) V^(x)| p efe < a fl (r) 

JU N \B r 

where 

a R {r)= f (p((\x\-2R-4:) + ) p dx = Lj N N [ t N ~ 1 V [[t - 2R - 4)+) p dt 

JR N \B r Jr 

and note that a R (r) — > for r — > oo. At the same time, let i(r) = sup[ r i then 
Vx € B r , |u n (a:) -tta'OOl < K r + 4 + 2-R)|u n (a;) - u n ,{x)\ 

so 

\v n {x)-v n i(x)\ p dx<u N r N l(r + 4 + 2R) p sup \u n {x) - u a \x)\ p < 

x£B r 

< uj N r N l(r + 4 + 2R) p d H (n, 0') P 

Summarizing, 

d{Q,n') p =f \v Q (x)-v Q ,(x)\ p dx+ [ \v Q (x)-v QI (x)\ p dx< 

JR N \B r J B r 

< a R {r) + oj N r N l{r + 4 + 2R) p d H {tt, 

Let eventually 

g R (s)= inf [a R (r)+u N r N l(r + 4 + 2R) p s] 

r>2R+i 



and note that it is concave and that lim^o g R (s) — 0; and let f R (s) — fy / g R (s p ). □ 
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Combining the two lemmas [3~9l and 1X71 we obtain that 

Theorem 3.10. The topology induced by d over the space X is equivalent to the one 
induced by dn- 

This implies that all properties of the Hausdorff distance listed in proposition ^. 181 
are valid for the distance d as well. 

3.1 Completeness and compactness 

By prop. 13.101 we know that (X, d) is locally compact. 
We now prove that it is complete: 

Proposition 3.11 (Completness). The space (X, d) is complete. 

Proof. Let f2„ be a Cauchy sequence; this means that, {vQ n }„ C N c is a Cauchy 
sequence: since LP is complete, VQ n — ► g in LP . It is well known (see e.g. thm IV.9 
in ID) that, up to subsequence that we indicate with {vk} k , there is also convergence 
Vk(x) — * g{%) for almost all x; let Uk(x)=ip~ 1 Vk(x) and u = tp~ 1 g; then Uk(x) — > 
u(x) on a dense subset, so by the lemma l2~T"5l u = uq where Q= {u — 0}. □ 

Summarizing, this and l3.10l implv that N c is a complete (that is, closed) and locally 
compact subset of LP . 

Remark 3.12. The above implies an interesting property of the subset N c of LP: it 
admits a small neighbourhood U on LP such that, for / € U, there is at least a v 6 N c 
providing the minimum of the distance inf v g^ e ||/ — u||. As far as we know, this 
minimum may fail to be unique. 

3.2 Shape analysis 

The family of distances is suitable for Shape Analysis: we can indeed prove 

Proposition 3.13. Let G = O(N) k R n be the Euclidean group of rotation and trans- 
lation; as in ( II Oi l, we can define the quotient metric by 

d q {[A],[B])=mid{ 9 A,B). (31) 

Then the above infimum is a minimum; so d q ([A], [B]) > when [A] ^ [B]. 
Proof. Choose a minimizing sequence {g n = (Rn, T n )} n£M , that is 

inf d(gA,B)= lim d(g n A,B)= lim d(R n A + T n , B). 

g£G Tt — >oo n — >oo 

Then {T n } nels4 must be bounded; in fact, let us assume by contradiction that \T n \ — > 
oo, then by (F22l we would have that 

d(A,B) < \\v a \\lp + \\vb\\lp 

and by <|23) that 

lim d{R n A + T n ,B) = \\v A \\ LP + \\v b \\lp, 
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so {g n } is not a minimizing sequence. This contradiction is generated by the assump- 
tion that {T n } is unbounded; then the translation part of every minimizing sequence of 
OTb must be bounded. By compactness we have that there exists a limit transformation 
g = (R, T) G K such that g n — > g and by continuity of d(fA, B) with respect of 
/ G G, we have that d(c?A, B) = d 9 ([A], [B]). □ 



3.3 d 9 and geodesies 

In this section we restrict p G (1, oo). 

Unfortunately d = d p>(p is not path-metric: 

Proposition 3.14. Given any two A,B el with A =/= B 

• then there is at most one A G (0, 1) such that Xva + (1 — X)vb G N c 

• consequently, by thm. \2.1l\ we have that d(A, B) < d 9 (A, B). 

Proof. It is immediate to show that f\ = Xva + (1 — X)vb assumes the value cp(0) 
only on the intersection of the two sets A n B for any A G (0,1). Then f\ G N c 
implies that fx = vadb- Let x e (An B) c such that va(x) ^ vb{x). We have that 
f\ x {x) 7^ f\ 2 {x) for any Ai, A2 G (0, 1) and Ai 7^ A2. Then there is at most one 
A G (0, 1) such that f\(x) = vac\b{x). □ 

So, to prove that the metric d admits minimal geodesies, we have to study d 9 as 
well; to this end, we prove two results. 

Proposition 3.15. If 

<p'(\x\) G L P (R N ) (32) 
then the space (X, d Pjt p) is Lipschitz-arc connected. 

Proof. Indeed, let j(t) — tft be the path that rescales 51 to the singleton {0}; we prove 
that 7 is Lipschitz. 

It is not difficult to prove that the map (t, x) Utn(x) is jointly Lipschitz. Then 
Utn(x) is differentiable at almost all t, x, and fix such a t, x; note that 

um(x) = tun (j) 

(as in eqn. d30l>): hence, taking derivatives w.r.t. x we obtain 

Vu m {x) = Vuq (Jj-^j 

while taking derivatives w.r.t. t we obtain 

d t um(x) = u n - j{Vu n (^j ■ x) = ~ (u m (x) - (Vu m (x) ■ x)) . 

Suppose now that x G" tfl and let y G tVl be a minimum distance point from x: then 

; - y 



um(x) = \x - y\ , Vutn(x) 
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so 



d t utn(x) = i ( \x- y\ - ( — — ■ x 

t \ \ f - v\ 

1 / \ , x -yy\ 

x-yy) = -{- l r ■ 7) (33) 



i|a; - j/| |a;- y\ t 

so if C -B r we obtain that \dtUtn(x)\ < r. If instead x E tfl and Wtn(x) is differen 
tiable at x then Vittn(x) = and dtUtn(x) = 0. 
To conclude (cf. 12.7b we compute 



llilli, = / \<p'(uta(x))\P\d t u m (x)\ p dx< r P J y{u m {x))\Ux (34) 
and we argument as in 13. 21 By Rem. 1.1.3 in 0, we conclude that 7 is Lipschitz. □ 

Remark 3.16. Asking that tp satisfy both ( fT9] l and (l32l l is equivalent to asking that 
G W 1,p . By using the equality in d34l l and in (l33l . it is possible to show that, 
for most compact sets, the rescaling is a Lipschitz path if and only if S W 1 '?. 

When X is Lipschitz-arcwise connected, the induced metric d 9 = {d v , v ) 9 is a finite 
metric. 

We can prove an equiboundedness result for d 9 (that is stronger than 13. 6b 

Proposition 3.17. Fix a compact nonempty set VL, and r > 0; then there is a K 
compact large such that for any closed set fl' satisfying d 9 (Q, Q') < r, then CI' C K. 

Proof. Let b(r) be defined in El Let d 9 (n,fl') < r, and 7 : [0,1] -> 7V C be a 
Lipschitz path (of constant L) connecting 7(0) = £1 to 7(1) = QI such that 

len 7 < d 9 (Cl,Cl') + 1 

up to reparametrization, we also assume that L < r + 2. Let n be large, so that 
(r + 2)/n < 6(r), and let K = fi + D rn (note that n only depends on r). Let 
= l{i/ n ) for i = 0, . . . , n; we know that 

d(Ai,A i+1 ) < d 9 {A ll A l+l ) < L/n<(r + 2)/n<b(r) 

since 7 is L-Lipschitz; so we apply recursively the proposition 13.61 on each Af. we 
obtain that, 

A l+1 c Ai + D r 

hence tt' C fl + D rn = K. □ 

The above results have many interesting consequences: 
Theorem 3.18. if 

<p'(\x\)eL*> 

then for any p > 0, 

D 9 (A,p)'M{A\d 9 (A,B) < p} 
is compact in the (I, d) topology; so 

• we obtain by Prp. \2.2\ that minimal geodesies do exist; 

• and b\ \2.4\ that the Geodesic Distance Based Averaging 



A = argmin A ^ d 9 {A, A,) 2 (35) 
of any given collection A\, . . . A n exists. 
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3.4 Variational description of geodesies 

In this section we restrict p 6 (l,oo). If j(t) is a Lipschitz path in N c , then it is 
associated to a function f(t, x) = v 7 m(x). 

Proposition 3.19. Suppose that t h-> /(i, -) is a Lipschitz path from t € [0,1] to 
L P (R N ); then, by \2.9\ for almost all t, f admits strong derivative -Jh in L P (R N ) (as 
was defined in eqn. (II 111 ). Moreover 

• / admits weak partial derivative dtf, and dtf = -i for almost all t. 

• If f admits a pointwise partial derivative hfor almost all t, x, then dtf = h. 

Proof. We extend f(t, x) = /(l, x) for t > 1, and f(t, x) = /(0, x) for t < 0; note 
that the extended f(t, •) is still Lipschitz in L P (R N ); then we define 

/. Ndcr/(* + t,x) - f(t,x) 

9r(t, X)= 

r 

so 

\\9r{t, x)||LP(RiV) < C 

where c is the Lipschitz constant of f(t, •); hence 

f \g T {t,x)\ p dxdt<c p 

This means that the family g T is bounded in L P ([Q, 1] x R N ), so we can find a sequence 
T n — > such that g Tii — > it) weakly, i.e. 

lim / / g Tn (t,x)ip(t,x) dxdt = / / w(t,x)ip(t,x) dxdt 
™ Jo Jk" Jo Jk n 

for all ^ £ C c °°([0, 1] x f). But 

' T (t,x)^(t,x)dxdt= [ f f(t, X ) ^ t ~ T,X l~^ t,X ^ dxdt 
JR N 

hence 



o Jr n Jo Jr n t 



lim / / g Tn {t, x)ip(t, x) dxdt = — \ \ f(t,x)dtif)(t,x) dxdt 
n Jo Jr n Jo Jr n 

by dominated convergence, so we conclude that / admits weak derivative, and the 
derivative is w. The relationship ( fT~3T > in L P (R N ), that is 

f(b,-)-f(a,-)= [ b f t dt 



implies that 



for all £ £ C c °°([0, 1]); but then setting ip(t, x) = £(t), we obtain that € = d t f. □ 
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This means that, for almost all t, we can represent the "abstract" derivative 2& by 
means of the weak derivative d t f{t, •) G LP(R N ). 

We use this result and eqn. ( TT4l > to express the length: 

len d 7 = / ||7(*)ll*= / Hf(t,x)\\ L pdt (36) 
Jo Jo 

So to find the minimal geodesic between two compact sets A, B, we need to minimize 
the above, with the constraint that /(0, •) = va, f(l, ■) = Vb, and, for any fixed t, 
(p^ 1 f(t, •) is a distance function. 

It is possible to prove (using a reparametrization lemma and Htielder inequality) 
that the geodesic is also the minimum of the action 

J(l)= [ \\d t f(t,x)\\ p LP dt= [ f \dtf(t,x)\*dxdt 
Jo Jo Jr n 

Equivalently, setting g(t, x) = u y f t -\(x), to find geodesies we can minimize 

J( 7 )= f I \^^(g))d t g(t 7 x)\Pdxdt 

Ja J\R N 

with the constraint that g(0, •) = ua, .9(1, ■) = Ub, and, for any fixed t, g(t, •) is a 
distance function. 



3.5 Tangent bundle 

Let p G (1, oo). We identify I with N c C L p , as by remark[3~4 
Given avE N c , let T V N C C LP be the contingent cone 



T V N C = \]imt n (v n — v)\tn> 0,v n G N c ,v n v\ = < AliiriT— - — H_ | A > 0, u„ 

L n - 1 I n |p n - V\\ LP 

where it is intended that the above limits are in the sense of strong convergence in LP . 

According to theorem lzTl if 7 : [a, b] — > N c is a Lipschitz curve then 7(£) G T 7 iV c 
for almost all t. 

In the following example we write explicitly the element of the contingent cone 
relative to a particular curve. 

Example 3.20. We fix £1 G I, and define the fattening f2 ( = £1 + Z?t for t > 0. We are 
interested in evaluating the derivative j(t). As previously done, we use the fact that 

uu t {%) = (uq(x) - t) + (37) 

and note that this map is jointly Lipschitz in (t, x): hence both un t (x) are va t (x) are 
almost everywhere differentiable. The pointwise derivative is given by: 

,. 1 r 1 f-<p'(u n (x) -t) for xg O t , 

w = lim - [v Qt+T - v n t \= i o (38) 

r^O T J 1 for X G fit ■ 

(note that the derivative may not exists for x £ dflt). If € LP then w G £ p , 

and it can be shown that 

to = lira - Ptffi t+T - un t l 

in the LP sense; then w is in the contingent cone. In particular, by Rem. 1.1.3 in J3], 
we obtain that the curve 7 is Lipschitz for t G [0, T]. 
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Unfortunately the contingent cone is not capable of expressing some shape motions 

Example 3.21. We consider the removing motion; to simplify the matter, let A be 
compact, and suppose that the origin is in the internal part of A; let A t =A \ B t be 
the removal of a small ball from A: then we can explicitly compute (for r > 0, s > 
small) 

\\VA r+s -vaAlp = 

+ 

< 

where L is the Lipschitz constant of (p(t) for small i, and we see that this motion is 
Lipschitz. If we try to compute 

\\vA t ~ V A \\ L p 

we notice that VA t — va = outside of B t : so the limit would be zero for x ^ 0. 

3.6 Riemannian metric 

Let now p = 2. The set N c may fail to be a smooth submanifold of L 2 ; yet we will, as 
much as possible, pretend that it is, in order to induce a sort of "Riemannian metric" 
on N c from the standard L 2 metric. 

We define the "Riemannian metric" on N c simply by 

(h,k)*(h,k) L * 
for h,k € T V N C and correspondingly a norm by 

|fc|V(M) 

Proposition 3.22. We will also argue that the distance induced by this "Riemannian 
metric" coincides with the geodesically induce distance d 9 . Indeed let 7 : [a, b] — » M 
be a Lipschitz curve in N c ; we may define the "Riemannian length " of the curve 

len R 7= J \~y\ds 

Then we define the "Riemannian distance " d R (x, y) as the infimum of\en. R 7 for all 7 
connecting x to y. But by eqn. d!4l >, len^ 7 — len 7 and d R = d 9 . 

3.7 Example: smooth convex sets 

We propose, as an example, an explicit computation of the Riemannian Metric. We 
fix p = 2, N — 2. Let O C R 2 be a convex set with smooth boundary; let y(8) : 
[0,L] — > (?f2 be a parametrization of the boundary, v(ff) the unit vector normal to dfl 
and pointing external to SI: then the following "polar" change of coordinates holds: 

^:R+x[0,L]^R\fi , i/>(p,0)=v(O)+pv(O) 



uj n N J t N - 1 (cp(0)-ip(s + r-t)ydt + 

lo n N J t^ 1 (ip(s - 1) - (p(s + r - t)Y dt < 
uj N r p L p (r + s) N 
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We suppose that y(9) moves on dfl in anticlockwise direction; so 

v=Jd s y , d as y = -KV , d s v = nd s y 

where J is the rotation matrix (of angle — 7r/2), k is the curvature, and d s y is the 
tangent vector (obtained by deriving y with respect to arc parameter). 

We can then express a generic integral through this change of coordinates as 

f( X )dx= f f f(lP(p,s))\l+ P K(s)\dpds 
R 2 \fl JR+ JdO. 

where s is arc parameter, and ds is integration in arc parameter. 

We want to study a smooth deformation of Q, that we call fl t ; then the border 
y(9, t) depends on a time parameter t. Suppose also that k(6) > 0, that is, that the set 
is strictly convex: then for small smooth deformations, the set fl t will still be strictly 
convex. By deriving 

d t d s y = d s (d t y) - d s y{d s y,d s (d t y)) = ■n v (d s (d t y)) 

where 

tt u (w)=w — v(y, w) 

is the projection of w parallel to v. Supposing now that p = p(t) as well, we can 
express the point ip(p, y) in a first order approximation as 



dip = ((d t y) + p'v + pjTT u (d s (d t y))jdt + (d e y + pd g v 



d6 



where moreover 



dgy + pdev) dQ = (d s y + pd s v) ds = (l + pn) d s yds . 

If y{9, t), p(t) are expressing a constant point x — ip(p, y), then dip — 0; we apply 
scalar products w.r.t. v and d s y to the above relations 

(u, (d t y)) + p' = , (d s y, (d t y)) - p(v, d s (d t y))dt + (1 + pn)ds = . 

Assuming that (dty) J- d s y, that is, (dty) = oiv with a = a(t, 9) £ R, we obtain 
the relationships 

, ds p(v,d s (av)) pd s a 

p = —a - — 



dt (1 + pn) (1 + pK) 
Now, for x <^ fl t , uq 4 (x) — p(t) hence 

h a =d t va t ( x ) = ~f' '(u nt (x))a 

whereas h a (x) = for x E tl t ; so h a is the vector in T V N C that is associated to a. 

Let us then fix two orthogonal smooth vector fields a(s)i , (s), 0(s)v(s), that repre- 
sent two possible deformations of dfl; those correspond to two vectors h a ,hp £ T V N C ; 
so the Riemannian Metric that we presented in Sec. I3.6l can be pulled back on dfl, to 
provide the metric 

(a, (3) = / h a (x)hp(x)dx = / h a (x)hf3(x)dx = 
Jr 2 Jr 2 \u 



on 



( l p'(p)) 2 (l + pn( S ))dp 



a(s)P(s)ds 
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that is, 



{a,j3)= / (a + b K (s))a{s)P(s)ds (39) 
Jan 



with 

Smooth sets If il is smooth but not convex, then the above formula holds up to the 
cutlocus. We define a function R(s) : [0, L] — + R + that spans the cutlocus, that is, 

Cut = {ip(R(s),s),se [0,L]} . 

ijj is a diffeomorphism between the sets 

{Go, s)s£ [0, L],0 < p < i?(s)} <-> R 2 \ (ft U Cut) 

moreover R(s) is Lipschitz (by results in 1 10|,| 12]). 
In this case the metric has the form 



(h,k) = 



on 



TS { V '{ P )f{l + P K{s))dp 







a(s)f3(s)ds 



4 Other Banach-like metrics of shapes 

The paradigm that we presented in the previous section may be exploited in other simi- 
lar ways; to conclude the paper, we shortly present some different embeddings (leaving 
to a future paper the detailed study of their properties). 

4.1 Signed distance based representation 

We may use the signed distance function bA, that was defined in (0, to define a metric 
of shapes: 

d'(A, B)=\\ip(b A ) - ¥>0b)||lp(rjv) 

in this case, we require that the function ip : R — > (0, oo) is monotonically decreasing 
and of class C 1 , and such that 

cp(\x\ - t) e ^(R^) Vt. (40) 

The resulting metric is slightly stronger than the one we studied in the preceding 
sections; in particular, 

Remark 4.1. Let T be the class of all finite subsets of R N ; this class is dense in X when 
we use the metric d P)(p , or the Hausdorff metric; but it is not dense when we use the 
metric d'. 

4.2 W 1 > p metrics 

Another interesting choice of metric is obtained by embedding the representation in 
W l p , forpg (l,oo) 
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We require that ip : [0, oo) — + (0, oo) be Lipschitz, C 1 and monotonically decreas- 
ing, and <^(|a;|) € W 1,P (R N ); for the case p < oo we are equivalently asking that 

t N -\ip(t)P + \ip'(t)\ p )dt <oc 

and this implies that lim^oo ip(t) = = lirrit_, 00 <p'(t). 

We add one last hypothesis, we assume that there is a T > s.t. <p(t) is convex for 

t e [T,o6\. 

Proposition 4.2. For any A compactwe have va € W 1 ' P (R N ). 

Proof. We already know by E2] that v A E L P (R N ). 

By hypotheses above, va is Lipshitz; and then, for almost all x, Vva — cp'(ua)Vua', 
where |Vu^| = 1 for almost all x £ A, while Vua = for almost all x £ A. We also 
know that when t > T, <p'(t) < 0, tp' is increasing and ip'(t) | 0. 

Let R > be large so that A C B R , then 

ua{x) >\x\—R 
and then when \x\ > R + T we obtain that 

y'{u A {x)) > <p\\x\-R) 

that is 

f \lp'{ua{x))\ p dx < [ \tp' {\x\ - R)\p dx < oo . 

Jr n \b r+t Jr n \b r+t 

At the same time, since va is Lipschitz, then J Br+t \ VvA\dx is finite. □ 
Definition 4.3. Given A,B el, we define 

di, p ,<p(A, B)=\\ip(u A ) - ip{u B )\\w 1 'P('R N ) 

We just state a simple property of this metric: 

Proposition 4.4. Let again T be the class of all finite subsets ofR N : this class is dense 
in X if and only ifip'(0) = 0. 

Indeed, fix A compact; let {x n } n be a dense subset of A; let A^= {xk \ k < n} be 
a finite subfamily; if(p'(0) — then Ak —>■ A according to d\ tPjV . 

If ip'(0) < it is easy to find examples where this does not hold: let N = 1, 
A = [0, 1], then \v' Ak (t)\ p dt —> |</(0)| p . 



Conclusions 

We have studied a metric space of shapes (1, d PtV> ); this space has a "weak distance", 
in that it has many compact sets, and geodesies do exist; but it can be associated in 
some cases to a smooth Riemannian metric, as we saw in eqn. d39t . Moreover, by the 
properties that we saw in sec. 12.21 (and in particular, by the properties of L p spaces 
for p e (1, oo) that we proved in Thm. !2.lTT i we can also hope that geodesies can be 
studied in the O.D.E. sense (altough possibly in a very weak sense). 

As we saw in the last chapter, the representation/embedding paradigm can be ex- 
ploited in many different fashions; we just conclude with one last remark. 
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Remark 4.5. The embedding of p o ua in W 2 ' p is not feasible: if A is smooth but is 
not convex, the second derivative of ua along the cutlocus is expressed by a measure 
(see 4.13 in 02)) and then tp o u A <£ W 2 > p . 
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